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ABSTRACT 

In order to understand the nature of the sources producing the recently uncovered 
CIB fluctuations, we study cross-correlations between the fluctuations in the source- 
subtracted Cosmic Infrared Background (CIB) from Spitzer/IRAC data and the unre- 
solved Cosmic X-ray Background (CXB) from deep Chandra observations. Our study 
uses data from the EGS/AEGIS field, where both datasets cover an ~ 8' x 45' region of 
the sky. Quantitatively, our measurement is the cross-power spectrum between the IR 
and X-ray data which we detect to be statistically significant and positive at angular 
scales ^ 20" where the source-subtracted CIB fluctuations in the Spitzer data are domi- 
nated by the clustering component. The cross-power signal between the IRAC maps at 
3.6/xm and 4.5^m and the Chandra [0.5-2] keV data has been detected with the overall 
significance of ~ 3.5(T and ~ 5a respectively. At the same time we find no evidence of 
significant cross-correlations at the harder Chandra bands. The cross-correlation signal 
is produced by individual IR sources with 3.6/im and 4.5/xm magnitudes rriAB^ 25-26 
and [0.5-2] keV X-ray fluxes ^ 7 x 10^"'^'' erg/cm^ /s . We determine that at least 
15 — 25% of the large scale power of CIB fluctuations is correlated with the spatial 
power spectrum of the X-ray fluctuations. If this correlation is attributed to emission 
from accretion processes at both IR and X-ray wavelengths, this implies a much higher 
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fraction of the accreting black holes than among the known populations. We discuss the 
various possible low- and high-z suspects for the discovered cross-power and show that 
neither local foregrounds, nor the known remaining normal galaxies and active galac- 
tic nuclei (AGN) can reproduce the measurements. These observational results are an 
important new constraint on theoretical modeling of the near-IR CIB fluctuations. 

Subject headings: 



1. Introduction 



Cosmic backgrounds contain emissions produced during the entire history o f the Universe in- 



cludi ng from objects individually inaccessible to telescopic studies (see review bv iKashlinskv et al 



20051 ). In different spectral regimes, the different cosmic backgrounds probe different sources ac- 
cording to their emission mechanisms. Thus cosmic X-ray background (CXB, ~[0.5-10] keV) probes 
either emissions by accreting back holes (BHs) or thermal X-ray emission from hot ionized gas, such 
as exists in galaxy clusters, whereas the cosmic infrared background (CIB) at the near-IR wave- 
lengths (1-5/im) is sensitive to stellar emissions. Correlations between structure in the IR and X-ray 
backgrounds could arise in two ways: they could be caused by one or more classes of sources that 
emit at both IR and X-ray wavelengths; or they could arise from separate classes of IR-emitting 
and X-ray-emitting sources that are found in association on large spatial scales. 

At the near-IR, the Galactic and Solar System foregrounds are substantial and, hence, must 



be kno wn to great accuracy when estimating the mean levels of the CIB. Thus iKashlinsky et al 
( 1996al lbl) and Kashlinskv &: Qdenwald (2000) pioneered the measurements of the CIB fluctuations. 



which circumvent many of the difficulties with the foreground subtraction. Indeed, the power spec- 
trum of the CIB fluctuations should reflect the clustering of the sources producing them. As the 
foreground galaxies get eliminated to fainter limits, the remaining source-subtracted CIB fluctua- 
tions would contain progressively larger fractions of the faint sources undetectable individually by 
current telescopic measurements. A particularly important class here are the sources associated 
with first stars epochs as the Universe gradually emerged from the "Dark Ages" . 

Current models predict the emergence of the first collapsed objects at redshifts z < 30 (see 
review by Bromm & Yoshida 2011). The expectation is that at these early times there will be an 
admixture of black hole (BH) populations which either formed by deaths of the first stars which 
formed in a top-heavy IMF, or formed by monolithic collapse of the primordial clouds. Although 
the first luminous objects and galaxies are too faint to observe on their own, fluctuations in the 
intensity of the cosmic infrared background (CIB) reflect the distribution of these early objects after 
foreground sources are removed to sufficiently faint levels. It has been suggest ed that these popula- 



tions may have left a meas urable signal in the mean CIB (ISantos 



20031 ) and its fluctuations (jKashlinsky et al.ll2004l : ICooray et al. II2004I ). There are intuitive reasons 



et a 



2002; Salvaterra &: Ferrara 



why CIB anisotropics from the early populations would be measurable: 1) first stars (and/or the 
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associated BHs) emitted a factor ~ 10^ more luminosity per unit mass than the present-day stellar 
populations, 2) their relative fluctuations would be larger because they span a relatively narrow 
time-span in the evolution of the Universe, and 3) they formed at the high peaks of the underlying 
density field which amplified their clustering properties. 

Intriguingly, there is now a substantial bod y of evidence suggesting that th e source-subtracted 
CIB fiuctuations, discov ered in recent Spitzer- (jKashlinskv et al.ll2005l . l2007al ) and AKARI-hased 
(|Matsumoto et aliboill ) studies, may arise from new populations which existed in the early Uni- 
verse. The residual CIB fiuctuations remain after removing galaxies to very faint levels and arise 
from populations with a significant clustering component, but only low levels of the shot noise 
(jKashlinskv et al.l l2007bl ) . This clustering signal ex ceeds, by a large and scale-dependent factor . 



the fluctuations produced by the remaining gala xies (jKashlinskv et al 



2005 



Helgason et al.ll2012l ). 



As suggested by jKashlinskv et al.l (j2005l . l2007bl ) thes e CIB fluctuations may or iginate in early pop- 
ulations. This found further support in a study by jKashlinskv et al.l (j2007d ) showing that there 
are no correlations between the source-subtracted IRAC maps and the faintest resolved sources 
observed with the HST ACS at optical wavelengths, which likely points to the high- 2; origin of the 
fluctuations, or at least to a very faint population not yet observed by other means. The high-z 
interpretation of the detected CIB anisotropics has received further confirmation in the recent Akari 
data analysis which measured source-subtracted CIB fluctuations to wavelengths as short as 2.4 fim 
and pointed out that th e colors of the fiuctuatio ns require their being produced by highly redshifted 
very luminous sources (jMatsumoto et al.|j201lj ). In an i mportant new step tow ard understanding 



the nature of these new populations, the r ecent study (jKashlinsky et al.jj2012l ) used Spitzer data 
from the SEDS program (jFazio et al.jj201ll ) and for the first time measured the source-subtracted 
CIB fiuctuations to ~ 1° showing that the amplitude of the CIB fiuctuations continues to diverge to 
more than 10 times that of known galaxies. The data indicate that these fiuctuations are produced 
by very faint sources and their angular spectrum is in agreement with the origin in early popula- 
tions spatially distributed according to the standard cosmological model at epochs coinciding with 
the first stars era. 

Such measurements, however, do not provide by themselves direct information whether the 
emissions in these new populations arise from stellar nucleosynthesis or BH accretion. If the sources 
producing these CIB fluctuations contained BHs in sufficient numbers, the latter sub-population 
would have contributed to the CIB fluctuations levels via accretion processes around the BHs. 
However, BH accretion produces a large fraction of emission in X-rays which could also produce 
a potentially identifiable component to the CXB. If the measured CIB fluctuations originate even 
partly from populations containing suflicient abundance of BH, the CXB component produced by 
them should correlate with the CIB providing a way to detect the BH population. 

Most of the CX B has now been resolved into points sources. Recent observations with Chandra 
(jLehmer et al.jj2012l ) resolved ~80-90% of the [0.5-7] keV energy band CXB into point sources. The 
majority of the sources contributing to CXB are AGN powered by accretion onto super massive 
black holes (SMBH). However, below the fluxes reached in deep Chandra observations, the most 
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numerous cla ss of sources are norma l galaxies whose X-ray emission is largely produced by X-ray 
binary stars. ICappelluti et alj (|2012l ) have shown, through angular fluctuation analysis, that 50% 
of the unresolved CXB is produce d by g alaxy groups and the remaining by galaxies and AGN. 
However, also in ICappelluti et alJ (j2012l ) it was suggested that if the large scale excess power 
observed in the GIB would be created by the accreting sources at z > 7.5, then up to 1/3 of the 
large scale GXB fluctuations could be produced by them without exceeding the observed power 
spectrum, and accounting only for a relatively small fraction of the total CXB flux. Indeed since 
high-z sources are expected to be highly biased they would produce backgrounds with a larger 
relative SF/F, and thus their fluctuation may be detectable despite a smaller contribution to the 
total GXB flux. 

Here we report the first direct evidence of the substantial X-ray er nission associated with 
the sources o f the GIB anisotropics uncovered in deep Spitzer /IRAC data (IKashlinsky et al.ll2005l . 
2007al . I2OI2I ) and briefly discuss the contributions to this signal from the various cosmological 
candidates. This result provides a major clue to the nature and epochs of the populations producing 
the source-subtracted GIB fluctuations. We detect correlations which indicate that at least 15 — 20% 
of the GIB is produced by objects with powerful X-ray emission. This proportion is much greater 
than among the known galaxy populations in the recent Universe. If the sources producing the 
sources-subtracted GIB signal are at high z, these findings suggest the possibility of substantial 
contribution from X-ray emissions to the reionization in the Universe. 

This paper is structured as follows: Sec. [2] discusses the data assembly for the EGS/AEGIS 
field observed by both Spitzer/IRAC and Chandra. Sec. [3] presents the results of the cross-power 
analysis identifying a highly statistically significant cross-power between the source-subtracted GIB 
and GXB. Finally, in Sec. H] we discuss the various possible low- and high- 2; contributors to the 
measurements. 



2. Data assembly 



2.1. X-ray data 



The primary X-ray data set used here is the deep Chandra AGIS-I AEGIS-XD survey (iGoulding et al 



2OI2I ) in the area overlapping with the SEDS IRAG survey in the EGS field. The relevant param- 
eters are listed in Table [U The field is located at Gelestial / Ecliptic / Galactic coordinates of 
(214.91°, 52.43°), (180.56°, 60.00°), (95.95°, 59.81°) and covers approximately 0.1 deg^. 

The Chandra X-ray Observat ory has a peak effectiv e area of 700 cm'^ at~1.2 keV and superb on- 
axis angular resolution of ~0.5" (iWeisskopf et al.ll2000l ). For imaging surveys, the X-ray telescope 
is generally coupled with a 16'xl6' GGD array, AGIS-I with an average energy resolution of ~130 
eV. The sensitivity window of Chandra covers the ^ [0.5 — 7] keV band, and since the GGD records 
the energy of the events it is possible to derive multi-band images with a single exposure. 
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The AEGIS-XD program consists in a series of 66 pointing in the central area of EGS field. 
For our purposes we used the deepest survey area and avoided regions near its East and Western 
border were the exposure drops dramatically. For every pointing we employed level-3 data produced 
for the Chandra -Source-Catalog with the most recent calibration database. The data have been 
cleaned of spurious events such as high ASCA grade cosmic rays as well as instrumental artifacts. 
Time intervals with high particle b ackground levels ha ve been removed. A detailed description of 



the data reduction can be found in lEvans et al.l (|201Cl ). Images have been created in the [0.5-2], 



[2-4.5] and [4.5-7] keV energy bands, respectively. The choice of this set of bands allows us to have 
the same number of counts (~1.3xl0^ cts ) and therefore the same statistical sampling in the three 
bands. In the same bands the exposure maps were computed at effective energies of 1.2, 3.2 and 
5.5 keV, respectively. Both images and exposure maps have been reprojected to match the IRAC 
maps at 1.2"/pix. Then, for each band, all the images and exposure maps have been combined 
into mosaic maps. The mean, cleaned exposure is 640 ks. Since we are interested in the source- 
subtracted CXB an important step in the data analysis is the point and extended sources removal. 
Thus, in order to remove as many sources as possible we performed a standard source detection 
in the [0.5-2] keV band and a combined [0.5-7] keV band by using the CIAO tool wavdetect with 
a threshold of 10^^, corresponding to < 5 spurious detections over the whole field of view. As a 
result we detected 252 and 303 sources down to fluxes of 7.0x10^^^ erg/cm^/s , and 1.1x10^^^ 
erg/cm-^/s , in the two bands, respectively. Note that the actual flux limit is model dependent. 
Here we assumed that the sources have a typical spectral index with r=2, and such a value can 
vary up to 20% if T changes of ±0.3. (No other sources than those detected in these two bands 
would have been detected in the [2-4.5] keV and [4.5-7] keV energy bands.) The actual CXB flux 
produced by detected sources is of the order 1.1 and 2.5x10^^ erg/cm^/s/sr in the [0.5-2] and 
[0.5-7] keV bands, respectively. These values carry an additional ±20% uncertainty because of the 
spectral model dependence. Our data are flux-limited in a position-dependent way and the values 
stated here are the average value of CXB resolved into point sources across the field of view. Our 
brightest source has a [0-5-2] keV fiux of the order of 5-6xl0~ ^^ erg/cm^/s , which is below the 



knee of the Log(N)-Log(S) distribution (jCappelluti et al.ll2009l ). thus a large fraction of CXB flux 



is not accounted as resolved flux mentioned above. Moreover since the X-ray maps used in this 
analysis are further masked for IR sources, the actual fraction of CXB resolved in our maps cannot 
be computed in a strayghtforward way. 

For every X-ray source we selected a circular region with the radius 5a of the PSF in o r der to 



allow an efficient source removal. Extended sources have been detected by lErfanianfar et al.l (j2012l ) 
by using a wavelet algorithm on scales of 32"-64" combined with optical red sequence identification. 
This procedure allowed us to mask clusters and groups of galaxies down to an estimated mass of 
~1O^^M0. The circular regions used here to mask extended sources enclose the projected r2oo 
radius, which ensures a highly efficient removal of the thermal X-ray photons contained in these 
groups. Since the Chandra PSF is energy dependent, in order to obtain a common mask for our 
X-ray maps, we selected for each the largest source extension a mong the three bands. The X- 



ray mask has been combined with the IR mask described below (jKashlinskv et al.l |2012| ) and is 
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depicted in the Fig(TJ The combination of the X-ray and IR mask left ~68% of the map pixels 
for fluctuation analysis via the FFT. The remaining counts are thus the CXB, plus the particle 
background recorded by the detector. The particle background has been subtracted by tailoring 
images taken by ACIS-I in stowed mode. Basically, ACIS was exposed when stowed outside the 
focal area. Since the particle background is not focused, the stowed image simply contains events 
due to particles. Such a background level however is not constant in time and thus one must find a 
recipe to renormalize the sto wed image to match the actual background level in the observations. 



Hickox fc: Markevitch I (j2006l ). showed that regardless of its amplitude, the particle background has 



a constant spectrum. In addition all the photons collected by Chandra in the [9.5-12] keV band 
have non astrophysica l origin (i.e. only particle events). Thus, the simple recipe proposed by 
Hickox &: Markevitch ] (&06) to compute the particle background level in each band is to scale 



the stowed images by the ratio Crfaja[9.5-12]/Csto«> [9-5-12], where Cdata and Cstow are total counts 
measured in the real images and in the stowed image, respectively. We have then subtracted the 
corresponding particle background image for each pointing. In addition, in order to compute the 
CXB fluctuation maps, we derived for every pointing and for every band, the mean CXB level map 
which is dependent on the off-axis angle because of vignetting. To do this, we created a map with 
a total number of counts equal to that of the real data outside the mask and distributing them 
according to the relative value of the exposure map. The count, mean-value and exposure maps 
have been then co-added in order to produce the final mosaicked Cx, {Cx) and E maps. The final 
fiuctuation image is then 6Fx = Cx/E — {Cx)/E. We also produced random noise maps drawn 
from two subsets of events. The events have been sorted in time and odd- and even-listed photons 
have been attributed to images a and b, respectively. These maps have the same exposure time and 
have been observed simultaneously so that effects of source variability are removed. In the same 
way as for real data we created a and b fiuctuation maps. The difference of these maps does not 
contain celestial signals or any stable instrumental effects. For this reason the |(a — 6) difference 
maps can be used to evaluate the random noise in the CXB fluctuations maps. Actually, the cosmic 
CXB fluctuation maps, 5Fx, have been produced by averaging the a and b data set, so that the 
auto- and cross-power spectra were evaluated on the ^(a + 6) maps. The fluctuation count-rate 
maps have been transformed into flux maps by applying the energy conversion factors (ecf ) listed in 
Table [U under the assumption that the average X-ray spectrum of the undetected sources could be 
represented by a power-law with r=2. Note that the actual spectrum of the sources contributing 
to the unresolved X-ray background is unknown since it is made by a blend of galaxies, clusters and 
AGN, and for this reason we have chosen the spe ctral model which reflects the mean spectrum of 



AGN and X-ray galaxies in the [0.5-7] keV band (|Ranalli et al.ll2005l : iGilli. Comastri. Sz Hasinger 
200?! ^. 



2.2. IRAC-based maps 



The Spitzer S pace Telescope is a 0.85 m diameter telescope launched into an earth-trailing 
solar orbit in 2003 (jWerner et al.ll2004l ). For nearly 6 years, as it was cooled by liquid He, its three 
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scientific instruments provided imaging and spectroscopy at wavelengths from 3.6 to 160 /im. In the 
time since the He supply was exhausted, Spitzer has continued to provide 3.6 and 4.5 fim imaging 
with its Infrared Array Camera(IRAC). The IRAC camera has a 5' x 5' field of view, and a pixel 
scale of 1.2", which slightly undersampled the instrument beam size of ~ 2" FWHM (jFazio et al 
20041 ). 



The procedure for map assembly is described in our previou s papers (IKashlinskv et al.ll2005l . 



2007al ) with an extensive summary including all the tests given in lArendt et al.l (|2010l ). Our IRAC 
mosaics are prepared from the basic calibrated data (BCD) product using the least-squares self- 
calibration procedure described by Fixsen, Moseley & Arendt (2000). The preparation and prop- 
erties of the IR data ob tained in the course of the SEDS program and used here are discussed in 
Kashlinsky et al.l (120121 ). The SEDS program was designed to provide deep imaging at 3.6 and 4.5 



Atm oy er a total area of about 1 square degree, distributed over 5 well-studied regions iFazio et al 



(j201ll ). The area covered is about ten times greater than previous Spitzer coverage at comparable 
depth. While the main use of the SEDS data sets will be the investigation of the individually 
detectable and countable galaxies, the remaining backgrounds in these data are well-suited for CIB 
studies, by virtue of their angular scale, sensitivity and observing strategies. Because of the suffi- 
ciently deep coverage with both X-ray and IRAC SEDS observation we have selected the Extended 
Groth Strip (EGS; Spitzer Program ID = 61042) field for this analysis. The field is located at 
moderate to high Galactic latitudes to minimize the number of foreground stars and the brightness 
of the emission from interstellar medium (cirrus). It also lies at relatively high ecliptic latitudes, 
which helps minimize the brightness and temporal change in the zodiacal light from interplanetary 
dust. The observations were carried out at three different epochs, spaced 6 months apart. At each 
wavelength, the frames are also processed in several different groups to provide multiple images 
that can be used to assess random and systematic errors. The noise is obtained by separating the 
full sequence of frames into the alternating even and odd frame numbers. Comparison of these "A" 
and "B" subsets, through construction of ^{A — B) difference maps, provides a good diagnostic of 
the random instrument noise because the A and B subsets only differ by a mean interval of ~ 100s. 

We also examined shallower (~ 3 hr integration) 5.8 and 8 /im observations of the EGS field 

Table 1. X-ray maps properties 



Band 


Nets" 




(Nph)/pix 




(CXBres ) 


ecf 










erg/cm^ /s 


xlO" *erg/cm^/s/sr 


xlO^i erg-l cm-^ 


0.5-2.0 keV 


233867 


133726 


0.23 


7xl0"i'^ 


1.10±0.08 


1.55 


2.0-4.5 keV 


216776 


137838 


0.23 






0.67 


4.5-7.0 keV 


201856 


134808 


0.23 






0.27 


0.5-7.0 keV 


652499 


406432 


0.69 


1.1x10-1'' 


2.5±0.19 


1.07 



"X-ray photon counts before masking. 
'^X-ray photon counts after masking. 
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Fig. 1. — From top to bottom : X-ray 0.5-2 keV fluctuation maps in counts rate units. The X- 
ray exposure map. IRAC-4.5/im fluctuation map. The IRAC 4.5//m exposure map. In all maps 
the black areas represent the mask. The EGS field is located at Celestial / Ecliptic / Galactic 
coordinates of (214.91°, 52.43°), (180.56°, 60.00°), (95.95°, 59.81°) and these subimages cover 
8'x45'. 

that were obtained during Spitzer's cryogenic mission (program ID = 8). However, even with 
application of the self-calibration, we find that the resulting images have background problems. 
Part of the problems are likely to be intrinsic and related to cirrus, i.e. thermal emission from 
interstellar dust. These observations cover a longer strip of the EGS than the SEDS observations. 
At the extreme end of the 8 /um image (~ 0.5° from the SEDS region) there is clearly diffuse 
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emission from cirrus, which is also evident in the IRAS 100 fim images and the LAB HI images 
(Neugebauer et al. 1984; Kalberla et al. 2005). The 5.8 //m data show additional background 
problems that are not correlated with the 8 iim data. These problems appear to be related to 
greater instability of the detector offset at 5.8 //m, which can be confused with temporal changes in 
the zodiacal light. Self-calibrating the 5.8 fim data without the subtraction of the estimated zodiacal 
light normally applied by the BCD pipeline provides a better, but still not satisfactory, result. The 
background issues at both 5.8 and 8 /zm may be compounded by the observing strategy. The 
SEDS strategy stepped across the full length of the field relatively quickly, and then accumulated 
depth by repeated observations, while the cryogenic observations accumulated the full depth of 
coverage at each pointing before moving on to another location along the EGS field. Because of 
these background issues and the higher noise levels in these data, cross correlations of 5.8 and 8 
fim emission with X-ray emission did not yield any significant results to present in this paper. 

The region selected for the joint CXB-CIB analysis is about ~ 8' x 45' in size. The common 
mask from IRAC 3.6/xm and 4.5/im and the X-ray was used, with about ~ 32% of the pixels lost 
to the analysis. An example of CIB fluctuation maps is shown in Fig. [TJ 

3. Fluctuation analysis 
3.1. Definitions 

The maps under study are clipped and masked of the resolved sources, yielding the fluctuation 
field, 5F{x). The Fourier transform, A(g) = J 6F{x) exp(—ix • q)Cpx is calculated using the FFT. 
The power spectrum in a single band n is Pnil) = (|A(^P)i with the average taken over all 
the independent Fourier elements which lie inside the radial interval [q,q + dq]. Since the flux 
is a real quantity, only one half of the Fourier plane is independent, so that at each q there are 
Ng/2 independent measurements of out of a full ring with Ng data. A typical rms flux 

fluctuation is y^q'^Pn{q)/2TT on the angular scale of wavelength 2n/q. The correlation function, 
C{6) = {6F{x) ■ 5F{x + 6)), is uniquely related to Pn{q) via Fourier transformation. If the fraction 
of masked pixels in the maps is too high, the large-scale map properties cannot be computed using 
the Fourier transform and instead the maps must be analyzed by direct calculation of C{6), which 
is immune to mask effects. In this study, the clipped pixels occupy ~ 32% of the maps which allows 
for a robust FFT analysis. 

We characterize the similarity of the fluctuations measured in different bands via the cross- 
power spectrum, which is the Fourier transform of the cross-correlation function Cmn(^) = {^Fm{x)- 
6Fn{x + 9)) . The cross-power spectrum is then given by Pmn{<l) = {^m{Q)^*n{l)) = T^miQ)T^n{Q) + 
Xm{q)'Fn{q) with TZ, X standing for the real, imaginary parts. Note the cross-power of real quantities, 
such as the flux fluctuation, is always real, but unlike the single (auto-) power spectrum the cross- 
power spectrum can be both positive and negative. 
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Fig. 2. — Left: Number of independent Fourier elements per bin that went into determining the power 
spectrum for each field. Right: The average number of photon counts per Fourier element adopted for 
determining the power spectrum as function of angular scale. 

The errors on the power have been computed by using the classical Poissonian estimators so 
that for the auto-power ap^ (q) = Pn{q) / -y/0.5 Nq and for the cross-power ap^^ {q) = Pm{q)Pn{(l) / Nq. 
These errors have been verified to be accurate to better than a few percent from comparison to the 
intrinsic standard deviation of the Fourier amplitudes at the various q. 



3.2. CXB povi^er spectra 

The analysis of the fluctuations of the CXB has been performed in the Chandra [0.5-2] keV, 
[2-4.5] keV and [4.5-7] keV bands. We evaluated the power spectra and their relative errors from the 
individual Chandra masked maps as well as from ^(a — 6) image. The final power spectrum of CXB 
fluctuations, Px-, is therefore evaluated as -Pi (^_,_^) — -Pi („_ with correspondingly propagated errors. 
The X-ray count maps, however, have an occupation number of <1 cts/pix, so the Gaussian behavior 
of their variance is not guaranteed especially at small scales. Correspondingly, we evaluated the 
mean number of photons per Nq/2 elements in the Fourier domain. The left panel of Figure [2] 
shows the number of independent Fourier elements, Nq/2, as a function of lii/q. The right panel 
shows the mean number of X-ray photons per element [i.e. Nets/ {Nq/2)] as function of angular 
scale, where Nets ~ 135000 (Table [1]). A limit of 20 cts/element is taken as a practical division 
between Gaussian and Poissonian regimes. The figure shows that below 10", X-ray counts are in 
the Poissonian regime, and therefore we limit our analysis of auto- and cross-power spectra to scales 
> 10" to avoid biases introduced by low-count statistics. 

In order to take into account the effects of sensitivity variation across the field of view, in every 
pixel the fluctuation field 6Fx{i) has been weighted by a factor E{i)/{E) where E{i) is the effective 
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Fig. 3. — Source-subtracted CXB fluctuation power spectra {top row) in the [0.5-2] keV, [2-4.5] keV 
and [4.5-7] keV energy bands and their relative cross-power spectra {bottom row). Open red circles 
and error bars represent the negative power points shown in absolute value for a better highlighting 
of the signal. 



exposure at the pixel i and {E) is the mean exposure in the field. The clipped and cleaned maps 
were Fourier transformed and power spectra evaluated. 

The binning of the power spectrum in angular scale is the same for all the energies sampled 
here. The relative sampling error (cosmic variance) on the determined power is [^Nq]~2, and so 
the power spectrum is not determined highly accurately at the largest angular scales (> 250") of 
the EGS field where lNg<10. The X -ray power spectra measured in the three X-ray energy bands 
are shown in Fig. [3j 

The [2-4.5] keV and [4.5-7] keV power spectra are noisier than that measured in the [0.5- 
2] keV band since in those bands the particle background is dominant with respect to CXB. In 
order to probe whether the source-subtracted maps at the different energy bands contain the same 
populations, we computed the cross-power spectrum between each pair of maps. This analysis shows 
that the cross-power spectra between the hard bands and the [0.5-2] keV band generally have lower 
amplitudes than the corresponding auto-power spectra, especially on smaller scales. This suggests 
that the population of sources producing the [0.5-2] keV CXB fluctuations is substantially different 
from that producing the hard X-ray CXB. 

It is important to emphasize, in the context of the discussion below (Sect. 4.), that the analysis 
of the CXB power and cross-power spectrum, carries an intrinsic source of uncertainty due to the 
contribution of the Galaxy. Thus the results shown in Fig. [3] for individual bands present an upper 
limit on the unresolved CXB fluctuations because they contain the contribution from the Galaxy 
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which is more prominent at the softest energies. Although our inspection of ROSAT Galaxy-CXB 
maps in this field does not show any well def ined structu r e, the actual shape of the Galaxy CXB 
power spectrum is unknown on these scales. Sliwa et al. ( 2001 ) measured the power spectrum of 
the ROSAT soft-CXB fluctuations and showed that its shape and amplitude is a strong function 
of the galactic coordinates. However, their measurements were obtained on scales larger than ~10' 
limiting any direct comparison to the CXB fluctuations in our field. Nevertheless, their Fig. 9 shows 
that the Galaxy component, at high galactic latitudes, is approximately white noise at sub-degree 
scales. Thus, while irrelevant for the CXB-CIB cross-power spectrum (see below), correcting for the 
Galaxy would reduce our estimate of the CXB auto-power particularly on the smallest scales. Since 
the Galaxy mostly emits below 1 keV, this could be the reason of a low-level of cross-correlation 
between [0.5-2] keV and [2-4.5] keV-[4.5-7] keV maps. In addition, even at the smaller scales, the 
CXB auto-power could suffer from an additional shot-noise component due to the wings of off-axis 
bright sources. Indeed at the largest off-axis angles the Chandra PSF is not precisely calibrated 
leading to less accurate masking of the brightest sources. Importantly, the reduced CXB fiuctuation 
power would make our conclusions about the unusually high contribution of accreting populations 
to the source-subtracted GIB sources even stronger. 



3.3. CIB power spectra 



In Fig. m we show the auto-power spectra of the IRAC 3.6//m and 4.5^m maps and their 
cross-power power spectrum. Th e CIB fluctuation spectr a evaluated in this work are in excellent 
agreement with those derived bv iKashlinskv et al.l (|2012l ) in the original EGS field even with the 
additional masking of X-ray detected sources. Power spectra with or without the additional X-ray 
masking agree within 5% on all scales, as shown by solid symbols and green lines. This is consistent 
with the populations responsible for the CIB fluctuation signal being unrelated to the remaining 
known galaxy or galaxy cluster populations in the field. 



3.4. CIB-CXB cross-povi^er spectra 

In order to establish if the fluctuations in the source-subtracted CXB and CIB maps have been 
totally or partly produced by a population of sources sharing the same environment (or even being 
the same sources), we performed the cross-power analysis and evaluated Pm,x{Q)- Since the X-ray 
and IR noise are uncorrelated, the cross-peer of the instrument noise contributions should average 
towards zero. 

The cross power-spectrum between IRAC 3.6/im and 4.5/im source-subtracted CIB fluctuations 
and Chandra [0.5-2] keV fluctuations are shown in Fig. [5l where we flnd a statistically significant 
cross-power. The same is plotted for IRAC 3.6fim and 4.5//m versus Chandra [2-4.5] keV and 
Chandra [4.5-7] keV, in Fig. [6l where we do not find statistically significant detection. 
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Fig. 4. — Left and central panels: The 3.6/xm and 4.5/xm CIB fluctuations power spectra, respec- 
tively in the EGS field. Right: 3.6//m vs 4.5/xm fluctuations cross-power spectrum in the EGS 
field. Green lines show CIB fluctu ations evaluated a f ter ap plying only the IR mask (instead of the 



IR-|- X-ray mask) to the data as in iKashlinskv et al.l (j201 



app 



We evaluated the overall significance of the cross-power by averaging the results on the whole 
angular range (10" < 27r/g < 1000"), computing the mean and its standard deviation. We also 
evaluated the significance from the actual dispersion of the un-binned data and found identical 
results. The mean power-spectra for every band pair investigated here are reported in Table EJ 
We find mean correlations at ~3.5cr and 5a significance between the [0.5-2] keV band and IRAC 
3.6 and 4.5 /xm bands respectively. The [2-4.5] keV and [4.5-7] keV bands do not show significant 
cross-correlation with IRAC band as shown in Tab. [2j We tested if the observed cross-correlation 
could have been produced by spurious instrumental features by cross-correlating the X-ray ^{a + b) 
maps with ^{A — B) IR maps and computed their average cross-power for every X-ray and IR band 
pair. The results of this analysis establish a noise floor for the cross-power spectra, and are plotted 
in Figs. [5] and [6] and listed in Table [21 We note that these X-ray signal vs. IR noise cross-power 
spectra are always consistent with zero and, as far as 3.6/im and 4.5^m versus [0.5-2] keV bands 
are concerned, the detected signal vs. signal cross-power is much larger than the noise floor. The 
noise floor confidence interval can be considered an estimate of the sensitivity of our analysis and 
a probe for systematic spurious power in the data. 

To further check the robustness of our results we also calculated cross-power spectra of our 
X-ray images with 1,000 random CIB fluctuation maps drawn by resampling the original maps. No 
statistically significant cross-power spectra were recorded. 
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Fig. 5. — (Top left) The fluctuations cross-power spectrum between IRAC 3.6/um 1/2(A-|-B) and 
Chandra [0.5-2] keV. (Top Right) the fluctuations cross-power spectrum between IRAC 3.6/im 
1/2(A-B) and Chandra [0.5-2] keV. The bottom row shows the same, but for IRAC 4.5/um and 
Chandra [0.5-2] keV. The labels in the plots list the average cross-power measured on the angular 
range 10"-1000" in units of 10"^" erg/cm^/s nW m~^ sr~^. Open-red-circles represent the negative 
power points and are shown on the positive quadrant for a better highlighting of the signal. 



4. Discussion 



4.1. Galactic, solar system, and instrumental foregrounds 



We begin by considering possible sources of spurious cross-correlation. In the IR bands, the 
most signi ficant foreground source of fluctuations would be the galactic cirrus. Yet, it was demon- 
strated bv lKashlinskv et al.l (j2012l ) that the bulk of the measured 3.6/im and 4.5/im power cannot be 
produced by cirrus. Cirrus emission is produced by dust in cold neutral and molecular clouds which 
cannot emit X-rays but can be effective absorbers of the so ft X-ray background leading to a neqatiy e 



contribution to the positive cross-power that is measured (jWang fc Yu 



1995 



Snowden et al.ll200d '). 



The Galactic X-ray emission from the hot phase of the ISM could play a role in the cross power, 
but there are several factors that limit its contribution to the cross power: 1) It is well known that 
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Fig. 6. — Same as Fig. [5] but for the [2-4.5] keV band (left panels) and [4.5-7] keV band (right 
panels). Note that in these hard X-ray bands the measured cross power is consistent with zero in 
all cases. 



the hot ISM mostly emits soft X-rays with energy <1 keV, and thus would be relatively weak even 
in the [0.5-2] keV band; 2) the Galactic X-ray background shows clustering on scales on the order 
of one degree, which is larger than the scales of interest here; and 3) the dust producing the IR 
emission is in the cold phases of the ISM, and thus should be anti-correlated with the hot ISM, 
leading to negative cross-power. 

Remaining faint Galactic stars c o uld, i n principle, contribute to the cross-power. However at 
high Galactic latitude iLehmer et al.l (|2012l ) showed that stars are a negligible component of the 
unresolved CXB. Moreover, the high level of isotropy of the GIB fluctuations works against the 
hypothesis of any Galactic sources as the possible sources of GIB fluctuations. 



Other possible sources of contamination discussed by iKashlinsky et al.l (j2012l ) are zodiacal 



Table 2. Mean Pir^x in units of 10 ^'^erg/cm^/s nW m ^ sr ^ computed over the [10"-1000"] 

angular range. 



Bands 
EcS 


0.5-2 keV 
1.2 keV 




2-4.5 keV 
3.2 keV 


4.5-7 keV 
2.3 keV 




{PlR,x) 


(PlA-B,x) 


{PiR.x 


> {P^A-B,x) 


{PiR.x) 


{PlA-B,x) 


3.6/xm 
4.5^m 


6. Oil. 7 
6.5±1.3 


-0.6±0.7 
-0.2±0.6 


3.2±2. 
2.2±2. 


6 0.1±1.4 
8 -0.9±1.4 


13.9±12.4 
9.9±9.4 


3.9±4.7 
0.0±4.4 



Note. — 



Bold text indicates the statistically significant results. 



- 16 - 



light and instrumental stray light whose contributions to CIB fluctuations were demonstrated to be 
negligible. At a low level, the IR zodiacal light may correlate with Solar System X-rays generated 
by solar wind charge exchange (SWCX). However, SWCX primarily produces very low surface 
brightness O VII emission around 0.54 keV where the Chandra effective area is very low and should 
not produce such a high signal. Moreover SWCX emission is time dependent and therefore since 
Spitzer and Chandra observed the field at different epochs, the signals are unlikely to show a 
correlation due to solar system effects. 

To summarize, our analysis points to an extragalactic origin of the positive cross-power spectra 
between the soft X-rays and the 3.6 and 4.5 fim background fluctuations. 



4.2. Extragalactic populations 

Several classes of extragalactic populations could contribute to the observed CXB-CIB cross- 
correlation. Below we briefly discuss the most obvious candidates for the emissions. More detailed 
interpretation will be worked out elsewhere, although it already appears that some of the candidates 
can be safely ruled out. 

It is important to emphasize that such populations lying at redshift z would have to contribute 
this signal in the X-ray energy range shifted upward by a factor {1 + z) in energy. Thus the 
populations beneath our detected signal would be emitting in the X-ray band of [0.5-2] (1 -|- z) keV 
at the effective rest energy of 1.2(1 -|- z) keV. 



As in iKashlinskv et al.l (|2012l ) we can express the common contribution of the sources to both 
IR and X-ray signals in terms of coherence, C{q) = j^^^(g)Px ■ The coherence can also be 
interpreted as the fraction of the emission due to the common populations so that C ~ CirCx^ 
where and Cx is the fraction of the emissions produced by the common population in the 
probed IR and X-ray bands. 

The source-subtracted CIB fluctuations are made up of two components: 1) small scales ($20" 
are dominated by the shot-noise from all sources (known and new) below the removal threshold, 
while 2) the larger angular scales refl ect CIB fluctuations produced by the clustering of the new 



populations iKashlinsky et al.l (l2007bl ). Thus the coherence between the two components of the 
fluctuations may be different depending on the different common levels of the populations producing 
the two terms. However, it is the larger scales, where the cross-power is due to clustering of the 
new populations common to both IR and X-ray emissions which are of greatest interest to interpret 
here. 

At large angular scales (iJ20"), where the clustering term dominates the CIB fluctuations 
spectrum, the coherence between 4.5/im and [0.5-2]keV is C ~ 0.02 — 0.05. The coherence at 3.6 
^m versus [0.5-2] keV is consistent with these values, although the cross-power is less statistically 
significant. 
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27T/q (orcsec) 27T/q (orcsec) 

Fig. 7. — (left): The IRAC 4.5/im fluctuation power spectrum (blue squares) compared with the 
4.5/um power from accreting sources (red circles), (right): The Chandra [0.5-2] keV fluctuation 
power spectrum (blue squares) compared with the power-spectrum of X-ray sources correlating 
with IR 4.5/im CIB. (red circles). The dotted line is the expected upper-limit for power spectrum 
for remaining X-ray normal galaxies with niA_B >25 at 4.5/im and X/0<-l. The dashed line is the 
expected upper-limit for power spectrum for remaining AGN with m^B >25 at 4.5/im and X/O<0. 
The continuous line is the sum of the AGN and galaxies fainter than m^^ ~25-25.5. 

Because the measured cross-power between the 4.5/im and X-ray data is highly positive, we 
plot in Fig. \7\ the CIB fluctuations produced by sources common to both the source-subtracted 
CIB at 4.5 fim and the [0.5-2]keV CXB, i.e. PciB.common = -PirxZ-^x = C x Pjr. This assumes 
that all of the CXB power spectrum is produced by these sources and implies a lower limit on the 
CIB fluctuations contributed by the common sources. The figure shows that 15 — 25% of the 
CIB power spectrum can be accounted for by these sources. The rightmost panel of Fig. [7] gives 
a similar plot for the minimal contribution, -PcxB, common ~ 

^iRx/-PiR = C X Px, to the CXB from 

the common populations. 

We must emphasize that the "common population" does not necessarily imply that the corre- 
sponding parts of the CIB and CXB are produced by the same physical sources emitting at both 
IR and X-rays. With the map resolution of a few arcsec we cannot resolve the individual point 
sources, especially if they are at high z. This is further amplified since the Gaussian regime of the 
X-ray maps is reached at angular scale of ~ 10" which subtend linear scale of ~ 0.1/i^^Mpc at 
z ~1 and this defines the scale of the individual "objects" in our analysis and the discussion below. 
Thus we cannot resolve whether the IR and X-ray emitters are one and the same or whether they 
are different sources that share the same environment at the relevant angular scales. Moreover, 
from the amplitude of the cross-correlation signal itself it is not possible to directly determine if 
the signal is produced by a single population of sources, or if it is produced by different popu- 
lations sharing the same environment. One can demonstrate through the Limber equation (see 
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e.g.. iKashlinskv et al.ll2005l : ICappelluti et al.ll2012l ) that the amphtude of the cross power spectrum 
depends on the product of the emissivity of the two population times the product of their large 
scale biases. This means that a population of highly biased sources correlated with a population 
of more weakly biased ones, in principle, could produce a substantial cross-power signal even if its 
volume emissivity (number density) is lower than that of the other population. 

For proper interpretation of the measured CIB-CXB correlation it is important to reiterate the 
limits imposed from the IR analysis itself. The so urces in the IRAC map s used here are removed 
down to the shot-noise PsN — 30 nJy- nW/m^/sr (IKashlinskv et al. 2012 ). which is equiva,lent t o 
sources removed to magnitudes of rriAB — 25 — 25.5 (jKashlinskv et al.ll2005l : iHelgason et al.ll2012l ). 
Therefore, for this discussion we adopt as the flux limit of ~ SOOnJy at the IRAC bands. 
Thus in order to account for the measured CIB fluctuation of 6F ~ 0.05 — 0.1 nW/m^/sr these 
sources must have projected angular number density n^5F / Sw^a x {SF/F)~^ where F ~ nS is the 
CIB level produced by them. The remaining CIB sources below the threshold would have to have 
(0.3 — 0.4)[(S'iini/300nJy)(5F/F)]~^arcsec~^ in order to explain the observed CIB at 3.6-4.5 ^m. 
Only the sources that can produce highly non-linear CIB fluctuations, 5F/F ^ 1 all the way to 
sub-degree scales, can have projected number density significantly lower than this. The CIB from 
such sources would, however, then exhibit a clear void-cluster CIB pattern contrary to what we see 
in the CIB maps. 

Consequently our measurements indicate that in order to explain the detected cross-correlation 
the sources producing them would have to account for ^ yfC or 15 — 25% o/ the CIB signal and 
be abundant enough to reproduce the required number density while accounting for the remaining 
CXB fluctuation. Correcting for the CXB fluctuations power from the Galaxy would reduce the 
CXB fluctuations and make our conclusions about the unusually high contribution of accreting 
populations to the source-subtracted CIB sources even stronger. Since X-ray sources are less 
abundant than optical/IR sources, a possible solution would lie in a numerous population of 3.6/xm- 
4.5/xm sources correlated with highly biased X-ray sources. 



4-2.1. Diffuse gas in clusters and WHIM 



As mentioned in the introduction, the r nain producers of the u nresolved soft X-ray CXB power 
are galaxy groups and the putative WHIM (ICen &: Ostrikerlll999l ). However the mass limit of our 
source detection allows us to exclude fron i our analysis ext ended sources with mass M > lO^^M©. 
Galaxy cluster scaling relations (see e.g., iPratt et al.ll2009l ) ensure in this case that these sources 
have a low kT (i.e. ^1 keV). Any sources correlated with clusters of galaxies at even marginally 
high z, would require the gas to be at temperature shifted upward by a factor of (1 + z), making 
the origin in this component even less likely, doubly so since the clusters/groups are expected 
to have colder gas at the early times. We therefore conclude that such a population cannot be 
responsible for these measurements. This is further confirmed by the fact illustrated in Fig. [4] 
which shows that the additional X-ray masking, which includes the resolved X-ray sources, does 
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not have any noticeable effect on the measured CIB fluctuations. Similarly, the WHIM, although 
it has never been significantly detected in emission, is expected to show a typical emission line 
dominated spectrum. Most of the emission is produced by H-, He-like O and Ne like ions, which 
emit at energies below < 700 eV. We also note that the diffuse sources producing the CXB peak 
at z ~ 0.1. If the observed cross-correlation arose at that redshift, then the IR sources would be 
a popu lation of still undetected numerous low luminosity (i.e. with L'^IO^ Lq, iKashlinsky et al 



2007bl ) galaxies, which further weakens this low-z hypothesis. 



4-2.2. X-ray emission in remaining "normal" galaxies 



X-ray binaries and supernova remnants are the main sources of X-rays from norm al galaxies. It 



has b een shown that this population emits X-ray with a typical spectrum n{E) oc E (jRanalli et al 



20051 ). therefore their emission could contribute to the whole energy range sampled here. Although 
it is not straightforward to determine the effective X-ray flux of galaxies with IR counterparts at 
mAB ~ 25 — 25.5, a useful tool to determine the effective X-ray brightness of galaxies below this 
magnitude limit is the X/0 ratio. In fact, it has been shown that for X-ray sources, the X-ray to 
optical/IR flux ratio assumes well deflned values according to the nature of the sources. 

The X/0 is defined as X/0 = log(/x//opt) = log(/x) + C + mw „„,/2.5. For the A.Bu m vs 
[0.5-2] keV band the constant C has a value of ~7.53 (extrapolated from lCivano et al. II2012I ). For 
observations with a depth comparable with ours, the value of X/O for normal galaxies is X/0 <C — 1 



( Xue et al. 



2011 



Thus X-ray galaxies with IR counterpart with mAB~25 — 25.5 should have the 
[0.5-2] keV flux / < 3 — 4 x lO'^^erg/crn ^/s , which is about one order of magnitude below the flux 
limit of the 4Ms CDFS (jXue et al.ll201ll ) and ~ 20 times fainter than our hmit for the EGS fleld. 



In order to determin e if these faint X-ray s ources could produce the observed fluctuations 
we adopted the recipe of ICappelluti et al.l (j2012l ) and computed the expected CXB fluctuations 
angular auto-power produced by the clustering component of these sources, i.e. the power in the 
X-ray bands from normal galaxies below z ~ 7.5 and /([0.5 — 2]keV) < 3 — 4 x lO^^^erg/cm^/s . 
This contribution, which is of the order of 2-3% of the total CXB power, is shown in Fig. [71 right 
and is systematically small compared to the measured power on scales [20 — 200]". 

Further problems with the proposed candidacy of the remaining normal galaxies are that 
1) the shot-noise component in the CIB fluctuations on small scales, which is dominated by the 
remaining normal galaxies, appears uncorrelated with that in the CXB as is shown by the drop in 
the correlated power at the smallest scales (see Fig. [7]), and 2) on large scales, which are dominated 
by the clustering component, the minimal C IB fluctuation s h own i n Fig. [7] appears larger than the 
normal gala ,xv component re c onstru cted by iHelgason et al.l (j2012l ) as displayed in the lower right 
of Fig. 9 of iKashlinskv et al.l (I2OI2I ). 



Thus normal galaxies could eventually be responsible only for a small part of the observed 
signal. 
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4-2.3. Remaining known AGNs 



Known AGNs may exhibit stron g IR emissions when they are surrounded by optically thick 
absorbing material (jElvis et alJll994l ). and/or if they are embedded in star forming regions. The 
contributions to the GIB from these sources at intermediate z ~ 2 — 4 would arise from the IR 
bump produced by hot dust with maximum temperature of ~ lO'^K. 

One should therefore consider whether known AGNs can be responsible for the observed cross- 
correlation between the source-subtracted GIB and CXB. A critical point in estimating their con- 
tribution to the measured cross-power is that the signal is produced by sources below the IR 
flux of Slim — 300 nJy at 3.6 and 4.5 //m which is fixed by the measured shot-noise level re - 



maining in the GIB maps (jKashlinsky et al 



2005 



Helgason et al. 



2012 



Kashlinskv et al.l 120121 '). 



Treister et al (2004, 2006) conducted a detailed Spitzer/IRAG-based census and modeling of the 
Type I and II AGNs in the GOODS region and their results show that one expects the total 
number density of Type I and II AGNs to be riACN — 6000deg~^ at the IR fluxes below Siim- 
This then results in the net GIB flux at 3.6 and 4.5 fiui from the remaining known AGNs to be 
/agn =i riAGN5'iim = 6 X lO'^MJy/sr or Fagn = 0.004 nW/m^/sr at 4.5 fim. Thus if the AGNs 
were to produce the me asured GIB signal of 5F ~ 0.05 nW/m^/sr at 4.5^m at sub-degree scales 
([Kashlinskv et al.ll2012l ). with their X-ray emissions accounting for the observed cross-power, the 
resultant GIB would have to have highly non-linear fluctuations on scales between 1' and 1° with 
5F/FZ 10. A possibility would be that the signal could be produced by a population of faint GIB 
galaxies correlating with highly biased high-z AGN. 



A new study IXue et al.l (|2012l ) reported a significant contribution to the unresolved GXB 
(~25%) at [6-8]keV from highly absorbed AGN with very faint optical counterpart (25 < m < 28 
at 0.85/^m). The quoted result is at 3.9a significance at [6-8] keV, while these populations are not 
detected below 4 keV. Thus they cannot be responsible for the observed effect since the correlated 
maps are all at energies effectively much below 6 keV. We further emphasize that only sources 
with 3.6 and 4.5^m fluxes below S ~ 300 nJy contrib ute to the meas ured fluctuations. Obscured 
AGN is the most abundant class among faint AGN (jHasingeii l2008l ) . They typically show very 
hard spectra and weak X-ray emission below ~ 3-5 keV. Since we did not detect a hard X-ray 
cross-power spectrum, these sources, if AGNs, would be either Type-I sources or high-z (z >2-4) 
obscured AGN, with their primary power-law component redshifted to the [0.5-2] keV band. 

Gappelluti et al (2012) calculated the expected clustering component of the angular auto-power 
spectrum produced by AGNs with IRAC 4.5Mm counterpart s with uiab >25-25.5. For the AGNs, 



the value of X/0 ~0 (jXue et al.ll201ll : iGivano et al. II2012I ). Therefore, in order to produce the 
observed cross -correlation they shoul d have [0.5-2] keV fluxes < 3 — 4 x lO^^^erg/cm^/s . By using 



the recipe of (jGappelluti et al.ll2012l ) we evaluated their expected angular auto-power under the 
assumption that they lie at z <7.5. Our prediction is shown in the right panel of Fig. [71 Its 
amplitude is of the order 7-8% of the total GXB fluctuations observed here. When added to the 
normal galaxies component this adds up to 10-11% of the total GXB fluctuations which is about 
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50% of the observed lower limit. 



4-2.4- New high-z populations 



Although no direct measurement of the redshift of the source-subtracted CIB fluctuations 
is yet available, there is now a significant body of evidence that the fluctuations likely originate 
at early times of the Universe's evolution: 1) The measured amplitude of the fluctuations can- 
not be accounted for by the low-luminosity e nd of the distribution of "ordinary" /known galaxies 



( Kashlinsky et al 



2005 



Helgason et alJl2012l ). 2) There are no correlations between the source- 



subtracted CIB maps at Spitzer wavelengths and HST/ ACS data out to 0.9 nm, which points 
to z > 7 — 8 for the populations producing the large scale excess signal unless the latter comes 
from new, and so far unobserved, very faint and more lo cal populations at AB mag ^i28 which 
have escaped the ACS detection (jKashlinskv et al.ll2007d ). 3) The pattern of the fluctuations is 
inconsistent with that of the g alaxy populations at rece nt times , and is consistent with the ACDM- 
distributed sources at high z (jKashlinsky et al.l 12007 bifl bold ). 4) The colors of the flu ctuations 
from 2 to 4.5 ^m are consistent with very hot sources at high-z (iMatsumoto et al.ll201ll ). 



If the X-ray signal comes from sources at high redshifts we clearly do not se e direct stellar 
emissions originating since massive metal-poor stars have T ~ (9 — 10) x lO^K ( Schaerer et al. 
20031 ). which is not hot enough to contribute to emissions in the observed X-ray bands extending 
to 7keV. Instead, if the signal originates from these sources, the contribution to the CXB signal 
would originate from thermal emission of the gas in accretion disks. 

We measure a coherence of C ~ 0.02 — 0.05. So if the BHs amongst the sources responsible 
for the measured source-subtracted CIB fluctuations produce the entire X-ray signal, they should 
account for ~ 15 — 25% of the signal produced at 4.5 /zm. If they contribute only a fraction of 
the X-ray fluctuations, their IR contribution would be even higher, but the measured cross-power 
suggests that 15 — 25%. At the lower limit of Cx the accreting sources (BHs?) would need to 
account for the entire CIB signal at 4.5 fim (i.e. CiR = !)• 

If the high-z sources are responsible for the detected cross-power, these early X-ray sources 
were present when the Universe was still partly neutral. Unlike UV photons. X-rays have the 
capability of multiple ionizations. If the sources responsible for the observed cross-correlation 
are at high-z sources, we are observing correlations between the visible (<4500A) and hard X- 
ray output of primordial accreting sources. Several authors suggested that early black hole X- 
ray feedback was necessary to reionize the Universe (jMadau et al.ll2004l : iRicotti &: Ostriker II2004I . 
20051 : ICiallongo et al.ll2012l ). and these results may suggest that the early universe was significantly 
radiated by hard X-rays which could have contributed to the reionization of the Universe. 
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4-2.5. New low-z sources 



If the CIB fluctuations that we have uncovered in Spitzer data arise from new populations 
at lower redshifts, say z ~ 2 — 4, they would have to originate in low mass (faint) system in 
order to account for the lack of correlations between Spitzer maps and ACS sources measured in 



Kashlinsky et alJ (l2007d ). The cross-power spectrum of CXB and CIB then requires that such a 



model would have to explain the existence of the significant BH emitters among these populations, 
sufficient to account for the observed contribution to the > [0.5-2] (1 + z) keV band emissions. 



5. Summary 

In the paper we presented the discovery of the statistically significant correlation between the 
3.6^m and 4.5//m source-subtracted CIB fiuctuations with the [0.5-2] keV CXB. Here we summarize 
the main results of the paper: 

• We detected a 3.5a to 5a significant cross correlation signal between the 3.6^m and 4.5^m 
source-subtracted CIB fluctuations and the Chandra-based [0.5-2] keV CXB fluctuations after 
masking X-ray detected sources and IRAC sources down to m^B 25-25.5. 

• With this dataset we do not find statistically-significant cross-power signal with the CXB at 
the harder X-ray Chandra bands ([2-4.5] keV and [4.5-7] keV). 

• The cross-power appears to be of extragalactic origin. 

• This result presents an important step in identifying the nature of the populations producing 
the source-subtracted CIB fluctuations discovered in Spitzer data. These populations must 
contain a signiflcant population of BHs which account for at least ~ 15 — 25% of the measured 
CIB signal. 
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